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Abstract: We present an analysis of the impact of the tree rings seen in the candidate sensors 
of the Large Synoptic Survey Telescope (LSST) on galaxy-shape measurements. The tree rings 
are a consequence of transverse electric fields caused by circularly symmetric impurity gradients 
in the silicon of the sensors. They effectively modify the pixel area and shift the photogenerated 
charge around, displacing the observed photon positions. The displacement distribution generates 
distortions that cause spurious shears correlated with the tree-rings patterns, potentially biasing 
cosmic shear measurements. In this paper we quantify the amplitude of the spurious shear caused 
by the tree rings on the LSST candidate sensors, and calculate its 2-point correlation function. We 
find that 2-point correlation function of the spurious shear on an area equivalent to the LSST field 
of view is order of about 10^^^, providing a negligible contribution to the 2-point correlation of 
the cosmic shear signal. Additional work is underway, and the final results and analysis will be 
published elsewhere (Okura et al. (2015), in prep.) 
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1. Introduction 

One of the most important studies in current cosmology is the constraint of the parameters of 
the standard cosmological model, including the equation of state of dark energy. Weak lensing 
is recognized as one of the most powerful methods to constrain the current cosmological model 
through the characterization of the mass distribution and evolution of the large scale structure (LSS) 
of the universe([l], [2], [3], [4], [5], [6]). 

As a consequence, several galaxy surveys that use the weak lensing of the LSS as one of their 
main techniques have started to operate (e.g. Dark Energy Survey^, HSC^, PanStarrs^, KIDS"*, etc.) 
or are being planned (e.g, LSST^, EUCLID^, WEIRST-AETA^). The cosmic shear analysis relies 
on the accurate measurement of shapes of millions or up to billions of background galaxies, to 
statistically estimate the small distortion (about 0.1%) caused by all the large scale structure in the 
universe. Therefore, it is crucial that all systematic errors comparable or larger than this threshold 
be corrected in order to accurately infer cosmological parameters. 

One of the instrumental features that has been identified in candidate CCDs for ESST and 
in detectors of other cameras (such as the Dark Energy Camera) is called tree rings ([7], [8], [9], 
[10]). Under uniform illumination, these appear as concentric circles whose center coincide with 
the center of the silicon wafer from which the detectors were cut. The origin of this effect lies in 

'http://www.darkenergysurvey.org/ 
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impurities in the silicon that generate spurious fields parallel to the surface of the CCD, transverse to 
the main electric field. The effective area of fhe pixels is modified and fhe phofogenerafed charges 
shiffed around, generating biases in phofomefric and asfromefric measuremenfs, if nof correcfed. 
In particular, fhe spafial variafion of fhe aslromefric displacemenfs disforfs fhe shapes of objecfs, 
producing a coherenf spurious signal fhaf mimics fhe coherenf disforfion produced by weak lensing. 
Therefore, if is imporfanf fo characferize fhe amplifude of fhe shear induced by fhe fhe free rings, 
fo defermine fheir impacf on cosmological paramefer esfimafions from weak lensing. 

This paper is organized as follows: in Section 2, we explain how we can esfimafe spurious 
shear from a concenfric displacemenf such as fhe one produced by fhe free ring effecl. In Secfion 3, 
we measure fhe spurious shear and ifs 2-poinl correlation by using fwo fypes of candidates CCDs 
for LSST CCDs. We summarize our findings and discuss furfher work in Secfion 4. 


2. Distortion by Concentric Displacements 


[7] (hereafter PBS 14) used images with uniform illumination (dome flat-field images, or flat fields) 
to characterize the tree rings of the 62 science CCDs of the Dark Energy Camera (DECam) [11] 
focal plane. After identifying the center of the tree rings in each CCD, radial profiles of their 
amplitudes as a function of distance from that center were created, assuming azimuthal symmetry. 
These profiles were directly used to quantify the photometric and astrometric displacements caused 
by the tree rings, and actual data from the science verification period of DES were used to constraint 
the astrometric displacements by the tree rings. The astrometric and photometric templates were 
incorporated into the photometric and astrometric (WCS) solutions to achieve residuals below the 
scientific requirements for DES. In this section we show how these kinds of displacements lead to 
distortions in the shapes of measured objects. 


2.1 Spurious Shear estimated from Concentric Displacements 

Eet us consider a displaced infinitesimal image between radius and and dd and 0^ + 50^ 

in polar coordinates, where the origin of the coordinate system is set to the center of the concentric 
displacement. The image is shifted from the original position and and 0^ and 0o + 8Qo 

by the tree ring displacement, which is written as a function d{rd). Eigure 1 shows the relation 
between the images. Then the edges of the image in the r direction relate to each other as follows: 

ro = rd-d{rd) ( 2 . 1 ) 

ro + 5ro = rd + drd-d {rd + drd). ( 2 . 2 ) 


The distortion can be obtained by taking the ratio of the length of the images in both the radial and 
tangential direction. Thus, the radial distortion is given by 

rp + 5ro - Vq _ ^ _ d{rd + 5rd)-d{rd) _ ^ _ ddjr) 
rd + 5rd - rd 5rd dr 

and the tangential distortion is given by 


, (2.3) 

r— 


ro5d _ ^ _ d{rd) 
rd5d rd 


(2.4) 
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Figure 1. Concentric displacement of an infinitesimal area by a distortion (i(r). 


For convenience, and without loss of generality, let us consider a situation where do = 0. Then we 
can arrange Eqs. 2.3 and 2.4 in matrix form, obtaining 

^ 0 \( 5rA^(l-K^^-yZ 0 U \ 

[roSeJ y 0 i-^JyrdSeJ y o + J yr^seJ 

(2.5) 

where and are the convergence and shear due to the tree-rings displacement, defined as 


^™(.) - 5 

rLd(f) = ^ 


dd{r) d{r) 
dr r 




(2.6) 

(2.7) 


Thus, the two components of the shear, y\ and 72 , at (r, Q) can be written as 

yl\r,e)+iyl\r,e) = yJo^Ar)e^‘^ 


(2.8) 


2.2 Spurious Shear estimated from Flux Modulation 

In practice, it is posible to directly measure the astrometric displacement d{r) caused by the tree 
rings by using sets of several dithered exposures of star fields in differenf phofomefric bands and 
consfrucfing “sfar flaf" images ([12], [13]). However, Ibis is nol an easy fask given fhaf fhe dis- 
placemenf is small (of fhe order of subpixels). PBS 14 demonsfrafed fhaf, fo firsl order, fhere is an 
relafionship befween fhe free-rings flux modulafion /(r) as measured by fhe flaf fields (which have 
higher S/N compared fo fhe one provided by fhe finife number of sfars in fhe sfar flafs) and fhe 
asfromefric displacemenf d{r), given by 

d(r) = -- [ dr'r'f{r'), (2.9) 

r 7o 
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The minus sign comes from the difference between the definitions of d{r). By using Eq. 2.9, we 
can rewrite Eqs. 2.6 and 2.7 as 

= -i/(r) (2.10) 

r^adir) = ( 2 . 11 ) 

where the last approximation is valid when the displacement is much smaller than the radius. The 
typical scales of the spatial displacements are of the order of subpixels, and the radii of the order 
of 100~1000 pixels, so the approximation is valid. This means that the convergence and shear due 
to the tree rings can be estimated as approximately half the value of the flux modulation. 

3. Spurious Shear and 2-point Correlation Function from tree rings on LSST CCDs 

In this section, we estimate the spurious shear caused by the tree-ring effect on two types of can¬ 
didate CCDs for the ESST project by using Eq. 2.11. It is important to keep in mind that this 
relation is true for an infinitesimal region, and therefore more realistic effects such as brightness 
distribution, pixelization, and PSE convolution need a more careful treatment. In particular, the 
PSE correction in weak lensing shear analysis should not be neglected when estimating the spuri¬ 
ous shear by the tree rings, because the tree-ring effect changes the ellipticity of the galaxies after 
PSE smearing. Therefore, in most cases, the PSE correction augments the spurious shear caused 
by the tree rings, e.g. in the situation of a circular PSE. However, in this paper we only estimate the 
impact on shear measurements due to the tree rings in a more simplified situation 

3.1 Flux modulation, displacement, and spurious shear by the tree rings 

We measured the tree-ring patterns of two different ESST candidate CCDs individually. Each one 
is fabricated by a different vendor, and we refer to them as “type 1" and “type 2" CCDs. We used 
flat-field images from a uniform light source of 750 nm wavelength. These flat images have a 4k 
X 4k pixel size and a standard deviation about 0.4%. 

Eirst we corrected the flat field images to eliminate other effects. We masked regions near the 
edge of the detectors and shadows due to dust particles in the optics. The light source in this test 
was not completely uniform, so we fit the slow variation with a 7th order polynomial. Other sensor 
effects are apparent in the image after correction due to pixel size variations and laser annealing. 
Though pixel size variation also will result in its own pattern of spurious shear, we removed both 
these effects by scaling the image. Though not ideal, this treatment is adequate for the present 
purpose. In future analysis we will study these effects in more detail. 

Next, we determine the center of the tree ring profiles in each CCD. We do this by selecting 
points on several rings and then fitting for their common center. The centers of the rings for type 1 
and 2 CCD lie at [4575,-375] pix [-167,4172] pix, respectively, where the origin of the coordinate 
system is at the bottom left corner of the CCD (see Eigure 4 for one of the CCD types). The position 
is near the corner, but slightly offset from it. 

The two panels in Eig. 2 show the measured tree-ring profiles of type 1 and 2 CCDs, respec¬ 
tively. Before measuring the profiles, the images were smoothed by a Gaussian function with a 
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5-pixel kernel. Sinee the seale of width of the tree rings peaks is of about 10 pixels, the smoothing 
should not reduee their amplitude by mueh. The typieal seale of the flux modulations is approxi¬ 
mately 0.01%, about 50 times smaller than that on the DECam CCDs (PBS14). The blue lines in 
Fig. 2 represent the tangential spurious shear ealeulated from the flux modulations by using Eq. 
2.11 (the spurious shear is about half value of flux modulation). The typieal amplitude seale of the 
spurious shear is of about 0.005%. 

Fig. 3 shows the astrometrie displaeement by the tree rings on both types of CCDs. The 
displaeement is ealeulated by using Eq. 2.9, whieh is an integral equation over the flux modulation. 
Therefore, large fluetuation errors from noise remain at large distanees. Sueh errors eould be 
eorreeted by using measured astrometrie displaeement from real stars ([7]). However, we ean say 
that the typieal amplitude is of about 0.001 pixels. 

3.2 2-point Correlation Function 

We now present the 2-point eorrelation of the spurious shear indueed by the tree-ring effeet on the 
ESST CCDs. First, we ealeulated the spurious shear pattern on eaeh individual CCD (Fig.4). The 
objeets(galaxies) we observe have finite size, so the spurious shear for them is not the same as the 
effeet on pixel seales. However, sinee the variation seale of the tree rings is larger than the objeet 
scale, we assume that the spurious shear effect does not change over the object scale.Therefore, the 
spurious shear for objects is almost the same as the effect on pixel scales, so we have calculated 
the 2-point correlation of the spurious shear for the latter. 

Fig. 5 shows the 2-point correlation function of the spurious shear on both CCD types. The 
typical value is about 10^ which is small, although at very short spatial scales we can see a 
slightly larger peak. This value is much smaller than the square of the typical scale of the spurious 
shear, (5 x 10^^)^ 10^^. The reason for the smallness is that the tree-ring oscillation gives 

images of objects an apparent spurious tangential and radial ellipticity in an alternate fashion, and 
therefore when calculating the correlation function the values cancel each other out. 

Fig. 6 shows the 2-point correlation of the spurious shear on type 1 and type 2 over an area 
equivalent to the full field of view of the ESST camera with a 100 pixel sampling scale. The typical 
scale for the amplitude decreases to 10^^^. In this calculation, we have assumed that all CCDs in 
the raft have the same tree ring profile but with a different orientation, since we have only tested 
one CCD of each type and four squares of silicon are cut out of the circular cross section of silicon, 
in such a way that the center of the boule, and of the tree ring pattern is at a corner of each square. 
Corners which have the center of tree-ring pattern are not same in all CCDs. 

4. Summary 

We have presented a study of the spurious shear induced by the astrometric displacement caused by 
the tree rings in ESST candidate CCDs. [7] (PBS 14) studied the tree rings on the DECam CCDs, 
and presented a model that estimates the astrometric displacement from the flux modulations by 
the tree rings as measured in flat-field images. Inferring the astrometric displacements in this way 
is relatively easier than direct measurements of the astrometric residuals from real data. 

We developed new formulae that calculate the spurious shear from the astrometric displace¬ 
ments. By using the fact that the tree-ring effect induces a concentric displacement, the spurious 
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Figure 2. 1-dimensional profile of the flux modulation(black) and the spurious shear(blue) caused by the 
tree rings effect on typel(top) and type 2(bottom) LSST CCDs. The tangential shear is estimated from the 
flux modulation by using Eq. 2.11. 



Figure 3. 1-dimensional prohle of the displacement caused by the tree-ring effect on typel (black) and 2 
(blue) LSST CCD. 


shear can be estimated from the derivative of the displacement function, where the displacement 
function is calculated from the flux modulation. We then measured the tree ring pattern on two 
LSST candidate CCDs, and calculated the induced shear. The typical value of flux modulation due 
to tree rings is 0.01%, about 50 times smaller than that for the DECam CCDs. In addition, the typ¬ 
ical amplitude of the 2-point correlation function in the field of view of the LSST is about 10^'^, 
which is small enough to be neglected in calculations that constrain cosmological parameters. The 
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Figure 4. Tangential spurious shear produced by the tree rings on typel LSST CCD, calculated from the 
flux modulation measured from the flat fields 
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Figure 5. Absolute values of the 2-point spurious shear correlations caused by the tree-ring effect on type 
l(top) and type 2 (bottom) LSST CCD. Black points mean parallel correlation and blue points mean cross 
correlation with 5 pixels sampling. The values oscillate around zero. 

correlation is much smaller than the square of the typical amplitude of the spurious shear due to 
cancelations of tangential and radial shears in the calculation. 

In this study we have measured only the impact of tree rings on the 2-point correlation function, 
but we plan to extend our analysis to include other sensor anomalies and the impacts of tree rings 
on cosmological parameters. 
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